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Abstract. We review properties of open chaotic mesoscopic systems with a finite 
Ehrenfest time te . The Ehrenfest time separates a short-time regime of the quantum 
dynamics, where wave packets closely follow the deterministic classical motion, from 
a long-time regime of fully- developed wave chaos. For a vanishing Ehrenfest time the 
quantum systems display a degree of universality which is well described by random- 
matrix theory. In the semiclassical limit, te becomes parametrically larger than the 
scattering time off the boundaries and the dwell time in the system. This results in the 
emergence of an increasing number of deterministic transport and escape modes, which 
induce strong deviations from random-matrix universality. We discuss these deviations 
for a variety of physical phenomena, including shot noise, conductance fluctuations, 
decay of quasibound states, and the mesoscopic proximity effect in Andreev billiards. 
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1. Introduction 

Technological advances of the past decade have made it possible to construct clean 
electronic devices of a linear size smaller than their elastic mean free path, but still 
much larger than the Fermi wavelength [1, 2, 3]. The motion of the electrons in these 
quantum dots is thus ballistic, and determined by an externally imposed confinement 
potential. On the classical level the dynamics can vary between two extremes, according 
to whether the confinement potential gives rise to integrable or chaotic motion [4, 5]. 
This classification is carried over to the quantum dynamics [6, 7, 8, 9], with most of the 
theoretical efforts focusing on the chaotic case. It has been conjectured [10, 11] that 
chaotic quantum dots fall into the same universality class as random disordered quantum 
dots. The latter exhibit a degree of universality which is well captured by random-matrix 
theory (RMT) [12, 13, 14]. Starting, e.g., from the scattering theory of transport [15, 16], 
RMT provides a statistical description, where the system's scattering matrix is assumed 
to be uniformly distributed over one of Dyson's circular ensembles [17, 18, 19]. This sole 
assumption allows to calculate transport quantities such as the average conductance, 
shot noise, and counting statistics, including coherent quantum corrections such as 
weak localization and universal conductance fluctuations [12, 13, 14]. RMT can also be 
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extended to the energy dependence of the scattering matrix by hnking it to a random 
effective Hamiltonian of Wigner's Gaussian universality classes [14, 19, 20]. This allows 
to investigate quantum dynamical properties such as time delays and escape rates [21], 
as well as the density of states of normal-superconducting hybrid structures [22]. 

RMT assumes that all cavity modes are well mixed, hence, that the system displays 
well-developed wave chaos. Impurities serve well for this purpose, since s-wave scattering 
couples into all directions with the same strength. Indeed, RMT enjoys a well-based 
mathematical foundation for disordered systems [20, 23, 24] . In baUistic chaotic systems, 
the scattering is off the smooth boundaries, which is not directly diffractive in contrast 
to the scattering from an impurity. On the other hand, classical chaos is associated 
to a strong sensitivity of the dynamics to uncertainties in the initial conditions, which 
are amplified exponentially with time ~ exp(At) (where A is the Lyapunov exponent). 
Hence, there are little doubts that RMT also gives a correct description for ballistic 
systems where this sensitivity yields to well-developed wave chaos. 

The precise conditions under which wave chaos emerges from classical chaos are 
however just being uncovered. Generically, RMT is bounded by the existence of finite 
time scales. In closed chaotic systems, for instance, spectral fiuctuations are known to 
deviate from RMT predictions for energies larger than the inverse period of the shortest 
periodic orbit [11]. For transport through open chaotic systems, classical ergodicity 
clearly has to be established faster than the lifetime of an electron in the system. 
Accordingly, the dynamics should not allow for too many short, nonergodic classical 
scattering trajectories going straight through the cavity, or hitting its boundary only 
very few times [25]. This requires that the inverse Lyapunov exponent and the 
typical time tb between bounces off the confinement are much smaller than the dwell 
time td, hence A~^, tb <S tq. In practice, these conditions are fulfilled when the openings 
are much smaller than the linear system size L. RMT universality also requires that 
A~^, Tb, and td are smaller than the Heisenberg time th = h/A (with A the mean 
level spacing). The condition Th ^ Tb guarantees that a large number of internal 
modes M = lnt[rH/TB] are mixed by the chaotic scattering (RMT is then independent 
of microscopic details of the ensemble [14]). The condition Th ^ r^, translates into 
a large total number of open scattering channels in all leads, A^tot = Int [th/td], such 
that details of the openings can be neglected. Together with the condition tb ^ td, 
this implies 1 <^ Atot <S M. The limit M — > oo, M/Atot — const is equivalent to 
the semiclassical limit of a small Fermi wave length Ap/-^ — > 0, all classical parameters 
being kept fixed. These requirements have been thoroughly investigated in the past 
[14, 24, 26]. 

More recently, following the seminal work of Aleiner and Larkin [27], it has become 
clear that a new time scale, associated to the quantum-to-classical correspondence of 
wave-packet dynamics (and in this sense, the validity of Ehrenfest's theorem), also 
restricts the validity of the RMT of ballistic transport. This time scale te, usually 
referred to as the Ehrenfest time, is roughly the time it takes for the chaotic classical 
dynamics to stretch an initially narrow wave packet, of Fermi wavelength Ap, to some 
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relevant classical length scale C (cf. Fig. 1). Since the stretching oc exp[At] is exponential 
in time, one has rg oc A~^ln[>C/AF] [28]. The Ehrcnfcst time poses a lower limit to the 
validity of RMT because wave chaos is associated to the splitting of wave packets into 
many partial waves, which then interfere randomly. In ballistic chaotic systems, the 
wave packet splitting is only established when initial quantum uncertainties blow up to 
the classical level. For shorter times, the quantum dynamics still bears the signatures 
of classical determinism, which is not captured by RMT. 

When Af is decreased, all classical parameters being kept constant - the very 
same semiclassical limit purportedly required for RMT universality - te becomes 
parametrically larger than tb and A~^, and indeed may start to compete with the dwell 
time td. One may thus wonder what is left of the RMT universality of open systems, 
and more generally of quantum effects in that limit. Indeed, there are many instances 
where quantum-to-classical correspondence at finite te leads to strong deviations from 
the universal RMT behavior. Such deviations are not only of fundamental interest, but 
also provide practical mechanisms to suppress or accentuate quantum properties. This 
short review provides a survey of the current knowledge of the quantum-to-classical 
correspondence in open ballistic systems, focusing on the deviations from RMT due to 
a finite Ehrenfest time. 

We start with a brief general classification of the Ehrenfest time for different 
physical situations such as transport, escape, and closed-system properties (section 
2). We then turn our attention to three specific applications where deviations from 
RMT universality occur once the relevant Ehrenfest time is no longer negligible. First 
(section 3), we discuss transport properties in a two-terminal geometry. Quantum-to- 
classical correspondence is refiected in the distribution of the transmission eigenvalues, 
and results in the suppression of electronic shot noise and the breakdown of universality 
for sample-to-sample conductance fiuctuations. Second (section 4), we discuss the decay 
modes (quasi-bound states) of the system. Escape routes faster than the Ehrenfest 
time give rise to highly localized, ballistically decaying quasi-bound states, while the 
density of long-lived quasi-bound states is renormalized according to a fractal Weyl law. 
Finally (section 5), we investigate the excitation spectrum of normal- metallic ballistic 
quantum dots coupled to an s-wave superconductor (the mesoscopic proximity effect). 
The presence of the superconducting terminal introduces a new dynamical process called 
Andreev reflection (charge-inverting retrorefiection) , which induces a finite quasiparticle 
lifetime and opens up a gap in the density of states around the Fermi energy. The size 
and shape of the gap show deviations from the RMT predictions when the Ehrenfest 
time is no longer negligible against the lifetime of the quasiparticle. Conclusions are 
presented in section 6. 

Because of the slow, logarithmic increase of Te oc InM with the effective size M of 
the Hilbcrt space, the ergodic semiclassical regime te > Td, Atd ^ 1 is unattainable 
by standard numerical methods. The numerical results reviewed in this paper are all 
obtained for a very efficient model system, the open kicked rotator [29, 30, 31, 32], which 
we briefiy describe in the Appendix. 
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Figure 1. Sketch of the dynamics of a wave pack:et in the phase space of an open 
system. The left and middle panels apply to the case of transport. The initial wave 
packet is maximally stretched along the stable manifold without substantial leakage 
out of the shaded rectangular area, which represents an opening of the system. After 
five bounces the wave packet returns to the opening, now being elongated along the 
unstable manifold (dashed lines; the sketches neglect the bending of the manifolds). 

(2) 

In the left picture, the transport Ehrenfest time is larger than the dwell time td- 
The returning wave packet still fits through the opening, with only minimal leakage. 
Hence, the particle leaves the system deterministically, as prescribed by the classical 
dynamics of the wave-packet center (dots). 

The middle picture corresponds to a more chaotic system (with larger Lyapunov 
exponent), resulting in < td. The stretching is stronger and the wave packet is not 
fully transmitted. In the subsequent dynamics, the partially reflected wave components 
will interfere randomly, which gives rise to wave chaos. 

In the escape problem, the initial wave packet can be squeezed more closely to the 
stable manifold, and the associated Ehrenfest time r^^^ is larger than the transport 
Ehrenfest time. This is illustrated in the right panel. 



2. Classification of Ehrenfest times 

The relevant Ehrenfest time depends on the physical situation at hand, but follows a 
very simple classification. Quantum-to-classical correspondence is maximized for wave 
packets that are initially elongated along the stable manifold of the classical dynamics, so 
that the dynamics first yields to compression, not to stretching (see Fig. 1). The initial 
extent of elongation along the stable manifolds is limited either by the linear width of 
the openings W, or the linear system size L, depending on whether the physical process 
requires injection into the system or not. In the same way, the final extent of the 
wave packet has to be compared to W or L depending on whether the physical process 
requires the particle to leave the system or not. For sufficiently ergodic chaotic dynamics 
(A~^, tb ^ Td, which implies W L) and in the absence of sharp geometrical features 
(besides the presence of the openings) , the resulting Ehrenfest time can be expressed by 
the three classical time scales and the Heisenberg time th [28, 33, 34]: 

P , . „-, I closed system, 

Tg"-* = In — I — I ; n = < 1 escape, (1) 
. B \ D/ J I 2 transport. 
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Here n gives the number of passages through the openings associated to the physical 
process. Expression (1) holds for two-dimensional systems such as lateral quantum dots 
[1, 2, 3] or microwave cavities [9], as well as for the stroboscopic one-dimensional model 
systems often used in the numerical simulations (tb is then the stroboscopic period; 
see the Appendix). Expression (1) also holds for three-dimensional systems (such as 
metallic grains) with two-dimensional openings when A is replaced by the sum of the 
two positive Lyapunov exponents. 

The difference between the three Ehrenfest times can be attributed to the additional 
splitting of a wave packet into partially transmitted and partially reflected waves at each 
encounter with an opening. Transport involves two passages via the openings. The first 
passage, at injection, determines the initial spread of the associated wave packet. The 
Ehrenfest time is then obtained by comparing the final spread to the width of the 
opening at the second passage, where the electron leaves the system. This results in the 
transport Ehrenfest time = ln[rHTB/T^]. The same Ehrenfest time also affects 
the excitation spectrum of normal-metallic cavities which are coupled by the openings 
to an s-wave superconductor, for which the relevant physical process is the consecutive 
Andreev reflection of the two quasi-particle species at the superconducting terminal (r^^ 
was actually first derived in that context [33]). In the escape problem, the electron is 
no longer required to originate from an opening. This lifts the restriction on the initial 
confinement of the wave packet and allows to squeeze it closer to the stable manifolds, 
as its elongation is not limited by the width of the opening but by the linear system size. 
Hence, the escape Ehrenfest time is larger than the transport Ehrenfest time by a factor 
td/t"b in the logarithm: r^^^ = A^^ ln[rH/TD]. This value is exactly in the middle of the 
transport Ehrenfest time r^^ and the conventional Ehrenfest time r^^ — A~^ln[rH/rB] 
for closed systems [28] , for which initial and final extents of the wave packet must be 
compared against the linear size of the system. 

The semiclassical limit is achieved for tr — > oo while the classical time scales 
A~^, Tb, and td are fixed. The Ehrenfest time then increases logarithmically with 
th/tb, which is oc L/\p for two dimensions and oc (L/Ap)^ for three dimensions. In 
this paper, we will denote this limit by M ^ oo while keeping M/Atot fixed, where 
M = Intfrn/ tb] is the effective number of internal modes mixed by the chaotic scattering, 
and Atot = Int[rH/TD] <^ M is the total number of open channels. 

3. Quantum-to-classical crossover in transport 

A rough classification distinguishes transport properties whose magnitude can be 
expressed by classical parameters from quantities that rely on quantum coherence [14]. 
Examples of the former class are the conductance G ~ pevp/V (where p is the electronic 
density in an energy interval eV around the Fermi energy, V is the voltage, and v-p is 
the Fermi velocity) and the electronic shot noise P ^ Pq — 2e^GV . The latter class 
is represented by the weak localization correction to the conductance and the universal 
conductance fluctuations, which in RMT are both of the order of a conductance quantum 
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Gq — e'^/h. In the presence of a finite Ehrenfest time, a much richer picture emerges: 
in particular the sample-to-sample conductance fluctuations are elevated to a classical 
level, while the shot noise is suppressed. 

The origin of these strong deviations from universal RMT behavior can be traced 
down to the distribution of transmission eigenvalues. We specifically consider transport 
through a chaotic cavity in a two-terminal geometry, and restrict ourselves to the case 
where the number of open channels leading to the electronic reservoirs are the same, 
A^L = -/Vr = N — A^tot/2. The scattering matrix 5 is a 2N x 2N matrix, written in 
terms oi N x N transmission (t and t') and reflection (r and r') matrices as 



r t' 
t r' 



(2) 



The system's conductance is given by G/Gq = TT{tH) = ^„T„ [15, 16], where the 
Tn G [0, 1] are the transmission eigenvalues of t'''t. In the hmit N ^ oo and within 
RMT, their probability distribution is given by [12, 13, 14] 

^™^<^' = ^VnT^ ■ ''^ 

Equation (3) requires that the standard conditions for RMT universality discussed 
in the introduction are met, A~^,tb <S td <S tr (hence, 1 <^ N <^ M). However, even 
when those conditions apply, it has recently been observed that strong deviations from 
Eq. (3) occur in the semiclassical limit M ^ oo [35]. This is illustrated in Fig. 2, which 
shows the results of a numerical investigation of the open kicked rotator (described in 
the Appendix). Instead of Eq. (3), the transmission eigenvalues appear to be distributed 
according to 

P^{T) = aPRMT(T) + [8{T) + 5{1 - T)] . (4) 

The presence of 5-peaks at T = and T = 1 in -Pq(T') becomes more evident once the 
integrated distribution I{T) = P{T')dT' is plotted. Prom Eq. (4) one has 

/„(T) = — arcsin Vt + ^-^(1 + 6i,t), (5) 

TT 2 

SO that /a(0) = (1 — a)/2 vanishes only for a = 1. For the data in Fig. 2, it turns 

(2) 

out that the parameter a is well approximated hj a ^ exp[— (r;^ + rB)]/rD] [35], 
with the transport Ehrenfest time given in Eq. (1). Hence, for a classically fixed 
configuration (i.e. considering an ensemble of systems with fixed A, tb, and td), the 
fraction / = 1 — a of deterministic transmission eigenvalues with T = 0, 1 approaches 
/ = 1 in the semiclassical limit M ^ oo, M/N — const. 

The emergence of classical determinism in the semiclassical limit reflects the fact 
that short trajectories are able to carry a wave packet in one piece through the system, 
provided that the wave packet is localized over a sufficiently small volume (see Fig. 1). 
Equation (4) moreover suggests that the spectrum of transmission eigenvalues is the sum 
of two independent contributions, precisely what would happen if the total electronic 
fluid of the system would split into two coexisting phases, a classical and a quantal 
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Figure 2. Left panel: Integrated probability distribution I{T) of transmission 
eigenvalues for open kicked rotators with td/tb = 25 and r^' {) (empty circles; 

('21 

distribution calculated over 729 different samples); td/tb — 5, and /tb — 0.16 (x; 

f2) f2) 

1681 samples), /tb = 1.5 (black diamonds; 729 samples), /tb = 2.8 (empty 

(2) 

squares; 16 samples), and Ytb = 4.1 (black triangles; 2 samples). The blue, violet, 
yellow, green and red lines give the distribution of Eq. (5), with a ~ 0.98, 0.81, 0.6, 
0.45 and 0.385 respectively. Right panel: Probability distribution P{T) of transmission 
eigenvalues for the same set of parameters as in the main panel (data for td /tb =5 and 
'T^ /tb = 1-5 have been removed for clarity). The blue solid line gives the universal 
distribution Prmt of Eq. (3) while the red dashed line corresponds to Eq. (4), with 
a = 0.39. Note that P{T) is symmetric around T = 0.5. 

one [36]. This splitting leads to a two-phase fluid model. It has been surmised that 
the quantal phase can be modelled by RMT [37], which results in an effective random- 
matrix model with renormalized matrix dimension aN. Since aN oc N^~^/^'^ oo 
in the semiclassical limit (see the related discussion of the fractal Weyl law in section 
4), effective RMT predicts that the universality of quantum interference such as weak 
localization and parametric conductance fluctuations is not affected by a flnite Ehrenfest 
time. This model is supported by a semiclassical theory based on the two-fluid model 
[38]. On the other hand, a stochastic quasiclassical theory which models mode-mixing 
by isotropic residual diffraction predicts that quantum interference effects are suppressed 
for a flnite Ehrenfest time [27, 39, 40, 41]. 

Numerical investigations on parametric conductance fluctuations [35, 40, 43] give 
support for the RMT universality of the quantal phase (see Sec. 3.2), and variants of the 
effective RMT model have been successfully utilized beyond transport applications (see 
Sees. 4, 5). On the other hand, while an earlier numerical investigation of the weak- 
localization correction [42] reported no clear dependence of the magnetoconductance 
5G = G{B = 0) — G{B = oo) on the Ehrenfest time, very recent investigations [39, 40] 
flnd a suppression of 6G for an increasing Ehrenfest time. The observed suppression 
is in agreement with the prediction SG oc exp[— /td] of a modifled quasiclassical 
theory [40] in which the suppression results from electrons with dwell time between 
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and . However, the quasiclassical theory cannot explain why the parametric 
conductance fluctuations arc not suppressed. It also does not yet deliver as many 
predictions beyond transport as the effective RMT (for quasiclassical predictions of 
the mesoscopic proximity effect, see Sec. 5). 

At present, both effective RMT as well as the quasiclassical theory have to be 
considered as phenomenological models, as they involve uncontrolled approximations. 
Clearly, a microscopic theory for the quantal phase which establishes the extent of 
its universality is highly desirable. This poses a considerable theoretical challenge 
considering that even in the limit of a vanishing Ehrenfest time a microscopic foundation 
for RMT in ballistic systems is only slowly emerging [44]. In this section, we focus on 
the consequences of the emergence of deterministic transport modes for the shot noise 
and the conductance fluctuations, as most of these consequences are largely independent 
of the precise degree of universality among the non-deterministic transport modes. 

3.1. Suppression of shot noise 

Shot noise is the non-thermal component of the electronic current fluctuations [45] 

/oo 
{51{t)5l{0)) exp{iujt) dt, (6) 
■oo 

where SI{t) = I(t) — I is the deviation of the current from the mean current /, and (. . .) 
denotes the expectation value. This noise arises because of stochastic uncertainties in 
the charge carrier dynamics, which can be caused by a random injection process, or may 
develop during the transport through the system. For completely uncorrelated charge 
carriers, the noise attains the Poissonian value Pq = 2e'^GV. Deviations from this value 
are a valuable indicator of correlations between the charge carriers. 

Phase coherence requires sufficiently low temperatures, at which Pauli blocking 
results in a regular injection and collection of the charge carriers from the bulk electrodes. 
The only source of shot noise is then the quantum uncertainty with which an electron 
is transmitted or reflected. This is expressed by the quantum probabilities < T„ < 1. 
In terms of these probabilities, the zero-frequency component of the shot noise is given 
by [46] 

Piuu = 0) = 2GoeVY,Tn{l - Tn). (7) 

This is always smaller than the Poisson value,which can be attributed to the Pauli 
blocking. 

In RMT, the universal value of shot noise in cavities with symmetric openings 
follows from Eq. (3), P{uj = 0) = ^Pq [12]. It was predicted by Agam et al. [47] that 
shot noise is further reduced below this value when the Ehrenfest time is flnite, 

p{oj = 0) = ^Po exp(-TE/TD). (8) 

The RMT value has been observed by Oberholzer et al. in shot-noise measurements on 
lateral quantum dots [48] . The same group later observed that the shot noise is reduced 
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Figure 3. Left and middle panels: dependence of the shot noise on the semiclassical 
parameter M ~ th/tb for two different Lyapunov exponents A and three different 
dwell times td of the open kicked rotator (as indicated in the plots). The solid line is 
the prediction from RMT. The dashed lines are obtained from a semiclassical estimate 
of the number of deterministic transport channels [37]. Right panel: rescaled data in 
a double- logarithmic plot, together with lines of slope I/Atd- Figures adapted from 
Ref. [30]. 



below the universal RMT result when the system is opened up (which reduces td/tb, 
not th/tb) [49, 50]. 

Equation (7) certifies that classically deterministic transport channels with Tn — 
orTn = 1 do not contribute to the shot noise [51]. Ref. [47] is based on the quasiclassical 
theory [27] which models mode-mixing by residual diffraction, and equates the Ehrenfest 
time with the closed-system Ehrenfest time rj^^ The discussion of the formation of the 
deterministic transport channels suggests that this has to be replaced by the transport 
Ehrenfest time [33]. Subsequent numerical investigations have tested this prediction 
for the open kicked rotator [30]. Results for various degrees of chaoticity (quantified by 
the Lyapunov exponent A) and dwell times td are shown in Fig. 3. The shot noise is 
clearly suppressed below the RMT value as the semiclassical parameter M is increased. 
The right panel shows a plot of — ln(4P/Po) as a function of In N'^/M ~ ln(THTB/4T^). 
The data is aligned along lines with slope 1/Atd. This confirms that the suppression of 
the shot noise is governed by , in agreement with the distribution (4) of transmission 
eigenvalues. 

3.2. From universal to quasiclassical conductance fluctuations 

Universal conductance fluctuations are arguably one of the most spectacular 
manifestations of quantum coherence in mesoscopic systems [52]. In metallic samples, 
the universality of the conductance fluctuations manifests itself in their magnitude, 
rmsG = 0{Go), independently on the sample's shape and size, its average conductance 
or the exact conflguration of the underlying impurity disorder. In ballistic chaotic 
systems, a similar behavior is observed, which is captured by RMT [12, 13, 14]. At the 
core of the universality lies the ergodic hypothesis that sample-to-sample fluctuations 
are equivalent to fluctuations induced by parametric variations (e.g. changing the Fermi 
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Figure 4. Variance cr^ (G) of the conductance vs. the rescaled effective Hilbert space 
size M/Mc in the open kicked rotator for parameters K e [9.65, 27.65], td /tb G [5, 25], 
and M e [128, 16384]. The scaling parameter Mc = 27r(TD/TB)^ exp(ArB) varies by a 
factor 70. The soHd and dashed hnes indicate the classical, sample-to-sample behavior 
oc M^, and the universal behavior (J^{G) = Gq/8, respectively. Inset: unsealed data for 
K = 9.65 and td/tb = 5 (circles), 7 (squares), 10 (diamonds), 15 (upward triangles) 
and 25 (downward triangles). Figure taken from Ref. [35]. 

energy or applying a magnetic field) within a given sample [52]. 

Three numerical works explored the quantum-to-classical crossover of conductance 
fluctuations [35, 40, 43]. Their findings are consistent with each other and support the 
conclusion that (i) the ergodic hypothesis breaks down once r^^^ is no longer negligible; 
(ii) under variation of a quantum parameter such as the energy, the conductance 
fluctuations stay at their universal value, independently on r^, /tq; and (iii) sample-to- 
sample fluctuations increase sharply above the universal value, cr^(G) oc GKM/Mc^, 
when M becomes larger than Mc — 27r(ro/rB)^ exp(ArB). 

Findings (i)-(iii) are illustrated in Fig. 4, which presents results obtained from the 
open kicked rotator model. They can be understood on the basis of the two-phase 
dynamical fluid discussed above. The deterministic transport channels are insensitive 
to the variation of a quantum parameter that influences the phase accumulated on an 
unchanged classical trajectories. However, once one changes the sample conflguration, 
all classical trajectories are scrambled and huge classical conductance fluctuations 
result. The size of the classical fluctuations is determined by the quantum mechanical 
resolution of classical phase space structures corresponding to the largest cluster of fully 
transmitted or reflected neighboring trajectories (see Ref. [37]), which yields the scahng 
with M (inset of Fig. 4) and (main panel of Fig. 4). When a quantum parameter is 
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varied, the conductance fluctuates only due to long, diffracting trajectories with t > , 
which build up the quantal phase. With the further assumption that the quantal phase 
is described by the effective RMT model, it follows that the parametric fluctuations are 

(2) 

universal, independently on (crosses in the main panel of Fig. 4). These conclusions 
are also supported by the observation in Ref. [35] that the energy conductance correlator 
F{s) = a^^{G) {6G{eo)6G{eo + e)) decays on the universal scale of the Thouless energy, 
oc h/Tj), independently on 



4. Decay of quaisi-bound states 

Suppose that the particle is not injected by one of the openings but is instead prepared 
(e.g., as an excitation) inside the system and then escapes through the openings (we will 
consider the case of a single opening with N^q^ = N channels) . Instead of the transport 
modes, this situation leads us to consider the decay modes of the system, determined 
by the stationary Schrodinger wave equation with outgoing boundary conditions. In 
contrast to the hermitian eigenvalue problem for a closed system, an open system with 
such boundary conditions features a non-self adjoint Hamilton operator with complex 
energy eigenvalues and mutually non-orthogonal eigenmodes, called quasi-bound states 
[19, 21, 53]. The imaginary part of the complex energy E = E' — i^ of a quasibound 
state is associated to its escape rate V (hence, all eigenvalues lie in the lower half of the 
complex-energy plane). These energies coincide with the poles of the scattering matrix, 
which establishes a formal link between transport and escape. Since RMT encompasses 
also the energy dependence of the scattering matrix, it delivers precise predictions for 
the escape rates and wave functions of the quasi-bound states. Hence, we are again 
confronted with the issue to determine the range of applicability of these predictions in 
light of the signatures of classical determinism observed in the short-time dynamics up 
to the characteristic Ehrenfest time for the escape problem. 

The universal RMT prediction for the escape rates can be obtained via two routes. 
The standard route relates the scattering matrix S to an effective M x M-dimensional 
Hamiltonian matrix H representing the closed billiard, and M x iV-dimensional matrices 
W that couple it to the openings [19], 

S{E) = l-2TxiW^{E - H + i-KWW^)-^W. (9) 

The superscript "T" indicates the transpose of the matrix. The poles of the scattering 
matrix are then obtained as the eigenvalues of the non-hermitian matrix H — inVVW'^. 
Assuming that if is a random Gaussian matrix, one can obtain detailed predictions of 
the density of these eigenvalues for arbitrary coupling strength [21]. For 1 <^ N <^ M, 
the probability density of decay rates is then given by 

nr) = -^e(r-r^i), (lo) 

where G is the unit step function. 

The second route to the distribution of decay rates is particularly adaptable for 
the case of ballistic dynamics. It starts from the formulation of the scattering matrix in 




Figure 5. Panel (a): Classical regions of escape after one to four bounces (color 
code) in an open kicked rotator with td/tb = 5 and K = 7.5 (Atb ~ 1.3). Panels (b) 
and (c): Husimi representations of two anomalously decaying quasi-bound states for 
M = 160. Panel (d): Husimi representations of a slowly-decaying quasi-bound state. 
Figure adapted from Ref. [34]. 



terms of an internal M x M-dimensional scattering matrix U{E), which describes the 
return amphtude to the confinement of the system [54, 55, 56]. For balhstic openings 
one has 

S{E) = PU{E)[1 - (1 - P^P)U{E)]-^P^, (11) 

where P oc W such that P^P is an idempotent projector onto the leads. The poles 
of the scattering matrix are obtained as the solutions of the determinantal equation 
det [1 - (1 - P^P)U{E)] = 0. 

In the semiclassical limit M oo, the matrix U{E) carries an overall phase 
factor exp(zz/(£^)) with phase velocity du/dE = tb/Ji, equivalent to a level spacing 
A = H/Mtb = h/r-n. In RMT, the distribution of the poles is obtained under the 
assumption that U{E) = Fexp{iETB/h) is proportional to a random unitary matrix F. 
The truncated unitary matrix (1 — P^P)F has eigenvalues /i = exp[—{iE' /h + 7/2)rB], 
where the decay constant 7 is related to the decay rate by F = (1 — e^'''^). The 
distribution of these decay constants is given by [57] 

^(7) = . ■ ,1 7^e(l - e-^- - tb/td), (12) 

4rDsmh (7x3/2) 

which is equivalent to Eq. (10). 

RMT does not account for escape routes with a lifetime shorter than the Ehrenfest 
time. Recently, it has been found that these routes induce the formation of anomalously 
decaying quasi-bound states, with very large escape rate F [34]. The semiclassical 
support of the associated wave functions is concentrated in a small area of phase space, 
(their total number is much larger than according to Weyl's rule of one state per Planck 
cell). For illustration. Fig. 5 shows classical regions of escape after a few bounces in 
the open kicked rotator, along with two examples of anomalously decaying quasi-bound 
states, which are both localized in the same region of classical escape after one bounce. 
These states are contrasted with a slowly decaying state, which displays a random wave 
pattern. 
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In order to discuss these observations, let us assume that the particle is initially 
represented by a localized wave packet Xo (such as sketched in Fig. 1). The evolution 
of this wave packet from bounce to bounce with the confinement is given by 

Xm = [(1 - P^PME)rxo. (13) 

When the final wave packet fits well through the opening the decay is sudden, hence 
not exponential at all. For such sudden escape the wave packets generated by the 
dynamical evolution all arc associated to rather special eigenstates of the truncated 
operator (1 — P'^P)U{E): If the escape occurs after n bounces, then 

[{l-P^P)UiE)r-^Xm^O (14) 

(neglecting the exponentially suppressed leakage out of the opening area). This 
corresponds to a highly degenerate eigenvalue — 0, hence, F = oo. 

Obviously, the (algebraic) multiplicity of this eigenvalue is at least m. However, 
in the semiclassical construction there is only one true eigenstate associated to this 
eigenvalue, namely, Xn-i- This deficiency is a consequence of the non-normality of the 
truncated unitary operator, for which the existence of a complete set of eigenvectors 
is not guaranteed. The degeneracy of the states is lifted beyond the semiclassical 
approximation, due to leakage out of the opening area. Hence, in practice one finds 
a complete set of eigenstates associated to this escape route, but the states are almost 
identical and hence are all supported by the same area in phase space. 

As a consequence of the strong overlap of the anomalously decaying states, Weyl's 
rule of covering the support of the states by Planck cells (of size ~ ^/M) cannot 
be used to estimate their number. The m states Xn, however, are semiclassically 
orthogonal, and Eqs. (13) and (14) imply that they span the same eigenspace as the 
nearly degenerate eigenstates (they provide a Schur decomposition). The orthogonality 
of these states reinstates the applicability of Weyl's rule. When we further observe that 
the semiclassical construction requires a reliable quantum-to-classical correspondence 
of the wave packet dynamics (mre < t^^), one can estimate the relative fraction / of 
baUistically decaying states by the probabihty to escape faster than Under the 

assumption of well developed classical ergodicity (tq ^ ''b,A~"^), this probability is 
given by 

/=l-exp(-4'V^D), (15) 

with the escape Ehrenfest time given in Eq. (1) [34]. 

Equation (15) has been tested numerically in the open kicked rotator by sorting 
all decay factors |//n| = 6xp(— 7„tb/2), n — 1, 2, 3, ... , M, according to their magnitude 
(see Fig. 6). The data approximately collapses onto a single curve when the relative 
index n/M is rescaled with exp(— t^^^/td). For small decay rates, the scaling function 
follows closely the RMT curve [57] 

ni\^i\)^M[l~^{l-\^i\Y'] (|/.r>l -tb/td), (16) 
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Figure 6. Ordered decay factors = exp(— 7„rB/2) for an open kicked rotator 
witli td/tb = 5 and K — 7.5 (Atb ~ 1.3) as a function of the relative indices n/M 
[panel (a)] and n/M [panel (b)], for M = 2™ x 80, m = 0, 1, 2, . . . , 8. The sohd hne in 
panel (a) is the RMT result of Eq. (16) with M = M. The solid hne in panel (b) is the 
same result with M given by Eq. (17-18). For the dashed line, this effective dimension 
has been fitted to the data. Figure adapted from Ref. [34]. 



where the matrix dimension is rescaled to 

M = Mexp(-4'V^D). (17) 
This equation can be rewritten as 

M = M'^-^/^^{tj)/tbY/^^, (18) 

which is precisely of the form of a fractal Weyl law [58, 59]. More generally, the 
exponent of M in this law is related to the fractal dimension of the repeller in the 
system. Equation (18) applies under the conventional conditions for RMT universality 
(tb, ^ td ^ th), for which the low fractal dimensions can be approximated by 
d = 1 — I/Atd + 0((ArD)~^) [60]. For a discussion of d and the scaling function outside 
this universal regime see the article of Nonnenmacher and Zworski in the present issue 
[59]. 



5. Quantum-to-classical crossover in the mesoscopic proximity effect 

We finally consider the situation of a ballistic metallic cavity in contact with a 
conventional superconductor, a so-called Andreev billiard [61, 62]. Compared to the 
normal billiards considered so far, the presence of superconductivity induces a new 
dynamical process called Andreev reflection, that is, retroreflection accompanied by 
electron-hole conversion [63] . This process prevents individual low-energy quasiparticles 
from entering the superconductor. 

For chaotic billiards it has been found that an excitation gap is formed as a 
consequence of the Andreev reflection, in that the Density of States (DoS) in the cavity 
is suppressed at the Fermi level. The energetic scale of this gap is the ballistic Thouless 
energy E'y = h/2Tj), where td is the average time between two consecutive Andreev 
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reflections [22] . For simplicity we consider the case of a single superconducting terminal 
with A^tot = N open channels at the Fermi energy. 

5.1. Bohr-Sommerfeld quantization versus random-matrix theory 

In an ergodic cavity, all classical trajectories except a set of zero measure eventually 
collide with the superconducting interface. Andreev retroreflection is perfect at the 
Fermi energy, where the hole exactly retraces the electronic path. For a nonzero 
electronic excitation energy E > 0, electron-hole symmetry is broken, and consequently 
there is a mismatch between the incidence and the retroreflection angle. For the lowest 
excitation energies oc in the limit Atd ^ 1, this mismatch is a small parameter. 
Consequently, all electron- hole trajectories become periodic. In the semiclassical limit 
where both the perimeter L of the cavity and the width W of the contact to the 
superconductor are much larger than Ap (hence M,N ^ 1), the semiclassical Bohr- 
Sommerfeld quantization rule relates the mean DoS to the return probability P{t) to 
the superconductor [22] , 

/*oo '^-^ n 1 fc~ 

p{E)=N / dtP{t)J2S " ' - 

Jo __n 



^-^^ ■ 2^ t 



(19) 



m=0 

The shift by 1/2 inside the (5-function is due to two consecutive phase shifts of 7r/2 at 
each Andreev reflection. This ensures that no contributions with an energy smaller than 

-E'min(^) = 7^h/2t emerge from trajectories of duration t. Since a chaotic cavity has an 
exponential distribution of return times P{t) oc exp[— t/ro] [64], Eq. (19) predicts an 
exponential suppression of the DoS in the chaotic case [65], 

Ntb {2tt/Ett))^ cosh(27r/ErD) 
TT sinh^(27r/ErD) 
The DoS can also be calculated in the framework of RMT. The excitation spectrum 
is obtained in the scattering approach from the determinantal quantization condition 
[14] 

Det[l + S{E)S*{-E)] - 0. (21) 

By Eq. (9), the scattering matrix S is then related to a Hamiltonian matrix H. For 
low energies, Eq. (21) can be transformed into an eigenvalue equation for an effective 
Hamiltonian [66] 

Det[E - i/eff] = , i/eff = ( _^ppT _H* ) ■ (22) 

Assuming that if is a random matrix, it has been found that the excitation spectrum 
exhibits a hard gap with a ground-state energy -Brmt ~ 0.6 £'t [22]. At first glance, 
both the Bohr-Sommerfeld and the RMT approach are expected to apply for chaotic 
cavities with A^^, Tb ^ t^. The hard gap prediction of RMT has thus to be reconciled 
with the exponential suppression (20) from Bohr-Sommerfeld quantization. 
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A path toward the solution to this gap problem was suggested by Lodder and 
Nazarov [67], who argued that the Bohr-Sommerfeld quantization is vaUd only for return 

(2) 

times smaller than the relevant Ehrenfest time, which was later identified with [given 
in Eq. (1)] [33]. For /td ^ 1, it was predicted that the hard RMT gap opens up at 

(2) 

an energy ~ in the Bohr-Sommerfeld DoS. A mechanism for the opening of this 

gap was soon proposed by Adagideli and Beenakker [68]. Constructing a perturbation 
theory, they showed that diffraction effects at the contact with the superconductor 
become singular in the semiclassical limit, which results in the opening of a gap at the 
inverse Ehrenfest time. More recent analytical and numerical works confirm that the 
solution to the gap problem lies in the competition between the Ehrenfest time and 
dwell time scales [29, 33, 69, 70, 71, 73, 74]. 



5. 2. Ehrenfest suppression of the gap 

At present there are two theories for quantizing Andreev billiards in the deep 
semiclassical limit. The first one proceeds along the lines of the two-phase fluid 
model and the effective RMT discussed in Section 3, but extended to take the energy 
dependence of the scattering matrix into account [70, 73]. The system's scattering 
matrix is decomposed into two parts, 

So{E) = S,i{E)(BS^^{E), (23) 

(2) 

where the classical part Sci{E) of dimension M(l — exp[— /tjj]) is complemented by 

(2) 

a quantal part Scs{E) of dimension Mexp[— /td]. The excitation spectrum hence 
splits into classical contributions originating from scattering trajectories shorter than the 
Ehrenfest time, and quantum contributions supported by longer trajectories for which 
diffraction effects are important. An adiabatic quantization procedure allows to extract 
the classical part of the excitation spectrum, while diffraction effects are included in 
the theory via effective RMT, SqmiE) — e-xp[iET^^ /h]S-RMT, where -Srmt is a random 
matrix from the appropriate circular ensemble, while the factor exp[iET^^ /h] accounts 
for the delayed onset of random interference. 

The second theory, due to Vavilov and Larkin [33], is based on the quasiclassical 
theory of Ref. [27], which models the mode mixing in the long time limit by isotropic 

(2) 

residual diffraction, with the diffraction time set to . Standard techniques based on 
ensemble averaging can then be applied. In the limit -C rD,thc predictions of both 
theories for the gap value (given by the smallest excitation energy eo) differ by a factor 
of 2 [62], 

eo , ari^^ f 2^/5-4 effective RMT, 

— — = 1 —■ S r- 1 ,1 (24) 

Ermt 2rD 1 V 5 — 2 quasiclassical theory. 

(2) (2) 

In the other limit ^ td, both theories predict eo = nh/2Tj^ [62]. In the transient 

(2) 

region ~ td, the two theories are parametrically different. These discrepancies 
have motivated detailed numerical investigations (based on the Andreev kicked rotator 
described in the Appendix) which we now review [29, 71, 73]. 
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Figure 7. Main plot: Dependence of the mean Andreev gap on the system size M, 
for open kicked rotators with td/tb = 5 and K = 14. Averages have been calculated 
with 400 (for M = 512) to 40 (for M > 5 • 10^) different positions of the contacts to 
the superconductor. The error bars represent the root-mean-square of cq. The dashed 
line is the RMT prediction and the solid line is a linear fit to the data points. Inset: 
Dependence of the mean gap on td/tb for K = 14 and M = 524288. The dashed 
line is the RMT prediction and the solid curve is given by Eq. (24), with coefficients 
extracted from the linear fit in the main plot. Figure adapted from Ref. [29]. 



We first show in Fig. 7 the systematic reduction of the excitation gap observed 
upon increasing the ratio r^^ /tj). The data corresponds to fixed classical configurations 
(dwell time and Lyapunov exponent) with variation of the semiclassical parameter M. 
The main panel is a scmi-logarithmic plot of e^/E-Y as a function of M G [2^,2^^], for 
t'd/tb = 5 and K — 14, well in the fully chaotic regime (Atb ~ 1.95). The data has 
been fitted to 

^ = 1 - [in(ArVM) - a'] , (25) 

as implied by Eq. (24) (the parameter a' accounts for model-dependent subleading 
corrections to the Ehrenfest time). We find a — 0.59 and a' — 3.95. Once a and a' 
are extracted, one obtains a parameter-free prediction for the dependence of the gap 
on tt)Itb, which is shown as the solid line in the inset to Fig. 7. We conclude that Eq. 
(24) gives the correct parametric dependence of the Andreev gap for small Ehrenfest 
times. Within the numerical uncertainties, the value of a conforms with the prediction of 
effective RMT. Similar conclusions were drawn in Ref. [74] from numerical investigations 
of Sinai billiards. 



5.3. Quasiclassical fluctuations of the gap 

The distribution P(eo) of the Andreev gap has been calculated within RMT in Ref. [75]. 
It was shown to be a universal function of the rescaled parameter (eo — E-^wv) / i 
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where A^v = 0.0Q8N^/^A gives the mean level spacing right above the gap in terms of 
the bulk level spacing A. Similarly, the standard deviation of the distribution is given 
by (7(eo) = 1.27AAr. 

The universality of the gap distribution is violated when the Ehrenfcst time is 
finite. As in the case of the conductance, the samplc-to-sample gap fluctuations are 
then dominated by classical fluctuations. In a simple approximation, the effective RMT 
model gives a qualitative prediction for the gap value in the crossover from a small to a 
large Ehrenfest time [70], 

A more precise form of the gap function was derived in Ref. [62]. Sample-to-sample 
fluctuations can be incorporated into the effective RMT model when one replaces the 
dwell time in Eq. (26) by the mean dwell time of long trajectories, that is one makes 
the substitution 

/OO / POO 

^^^dttP{t) I |(,,dtP(i). (27) 



This was done in Ref. [71]. The result with the correct gap function from Ref. [62] is 
shown in Fig. 8 [72]. It is seen that the gap fluctuations are greatly enhanced to the 
same order of magnitude as the gap itself. It was indeed found that oie^) becomes a 
function of td only in the hmit of large M. Prom Fig. 8, correlations between sample- 
to-sample variations of eo and (i)^^ are evident, clearly establishing the classical origin 
of the sample-to-sample fluctuations in the large M regime. 

Together with the average value and sample-to-sample fluctuations of eo, additional 
numerical evidence for the validity of the two-phase fluid model in Andreev billiards was 
presented in Refs. [71, 73]. Most notably, the critical magnetic field at which the gap 
closes was found to be determined by the competition between two values, B'^^ and B^. 
These fields correspond, respectively, to the disappearance of the gap for the quantum, 
effective RMT part and the classical, adiabatically quantized part of the spectrum. 
Moreover, Ref. [73] showed how most of the full density of states at finite t^^V'^d can 
be obtained from the effective RMT model. It would be desirable to have similar 
predictions from the quasiclassical theory that could be checked against numerical 
data. The excellent agreement between the numerical data and the predictions from 
the effective RMT model in Andreev billiards only adds to the intriguing controversy 
about the universality of the quantal phase in the two-fiuid model, which we encountered 
at various places in this review. 



6. Summary and conclusions 

We gave an overview over recent theoretical and numerical investigations which address 
the emergence of quantum-to-classical correspondence in mesoscopic systems with 
a finite Ehrenfest time. By now there is overwhelming evidence that the quasi- 
deterministic short-time dynamics up to the Ehrenfest time jeopardizes the universality 
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Figure 8. Quantum mechanical gap values eo of the Andreev kicked rotator as a 
function of the position picad of the center of the interface with the superconductor, for 
parameter values M = 131072, td/tb = 5, K = 14. The solid line uses effective RMT 
to relate the gap fluctuations to the fluctuations of the mean dwell time {t) » of long 
classical trajectories, defined in Eq. (27). Figure courtesy of Marlies Goorden [72]. 



otherwise exhibited by quantized open chaotic systems. This was illustrated in the 
discussion of three different physical situations: transport, decay of quasi-bound states, 
and the mesoscopic proximity effect. 

While there is consensus over the role of deterministic transport and decay modes, 
there are two competing theoretical frameworks with different predictions for the degree 
of wave chaos in the long-time dynamics beyond the Ehrenfest time, namely, the effective 
random- matrix theory [37] and the stochastic quasiclassical theory [27]. Both theories 
incorporate the deterministic short-time wave-packet dynamics in similar ways, and 
correctly explain the suppression of shot-noise power as well as the emergence of classical 
sample-to-sample fluctuations in the semiclassical limit. However, the two theories 
model the long-time dynamics in different ways, namely, via a random-matrix of reduced 
dimension or via residual diffraction. Consequently, the effective RMT predicts that 
quantum interference corrections like the weak-localization correction and parametric 
conductance fluctuations stay universal deep in the semiclassical limit, while the quasi- 
classical theory predicts a suppression of these effects in this limit. Moreover, there is 
conflicting numerical evidence for these coherent effects, with all the indications that 
more surprises are likely to be uncovered. In view of this intriguing situation, more 
theoretical and numerical efforts to uncover the limits of universality in mesoscopic 
systems, while challenging, are clearly desirable. 
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Appendix: open kicked rotators 

The logarithmic increase of the Ehrenfest time with the effectuve Hilbert space size M 
requires an exponential increase in the latter to investigate the ergodic semiclassical 
regime te ^ td, Atd ^ 1, in which deviations from RMT universality emerge due to 
quantum-to-classical correspondence. The numerical results reviewed in this paper are 

all obtained for a particular class of systems, the open kicked rotator [29, 30, 31, 32], 
for which very efficient methods based on the fast-Fourier transform exist. Combined 
with the Lanczos exact diagonalization algorithm, as first suggested in Rcf. [76], these 
methods allowed to reach system size in excess of M = 10® for the Andreev billiard 
problem [29]. 

The classical dynamics of the closed system are given by a symmetrized version of 
the standard map on the torus x,p E [0, 27r] [4], 



Each iteration of the map (28) corresponds to one scattering time Tb off the boundaries 
of a fictitious cavity. The dynamics of this system ranges from fully integrable {K = 0) 
to well-developed chaos [K > 7, with Lyapunov exponent A ~ ^ ln(i^/2)]. 

The map (28) can be quantized by discretization of the space coordinates Xm — 
27rm/M, m — 1,...M. The quantum representation is then provided by a unitary 
M X M Floquet operator F [8] , which gives the time evolution for one iteration of the 
map. For our specific choice of the kicked rotator, the Floquet operator has matrix 
elements 



The spectrum exp{—iEnTB / h) of F defines a discrete set of M quasienergies En G 
[0, /t/r^) with an average level spacing A = Ii/Mtb- 

For the transport problem, the system is opened up by defining two ballistic 
openings via absorbing phase space strips [xl — Sx,xl + Sx] and [xr — dx,XR -\- Sx], 
each of a width 2Sx = tttb/td. Much in the same way as in the Hamiltonian case [19], a 




(28) 



F 



exp{-(iMX/47r)[cos(27rm/M) + cos(27rm7M)]} 
X ^exp[2ml{m - m')/M] ex.p[-{m/2M)f]. 



(29) 
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quasienergy- dependent 2A'" x 2A'" scattering matrix can be determined from the Floquet 
operator F as [56] 

S(E) = P[ex.p(-iETB/h) - F{1 - P'^P)]-^FP'^. (30) 

The 2A'" X M-dimensional matrix P describes the couphng to the leads, and is given by 

p ^fl if n = me {(a:R,L-5)M/27r,a:R,L + 5)M/27r} 
I otherwise. 

The number of channels in each opening is given by = Int [dM/ir]. An ensemble of 
samples with the same microscopic properties can be defined by varying the position of 
the two openings for fixed Td/tb and or by varying the energy E. 

For the escape problem, P couples only to a single opening, and the quasibound 
states are obtained by diagonalization of the truncated quantum map (1 — P^P)F. 

So far we described particle excitations in a normal metal. In order to model 
an Andreev billiard [29], we also need hole excitations. A particle excitation with 
energy E^n (measured relatively to the Fermi level) is identical to a hole excitation with 
energy —Em which propagates backwards in time. This means that hole excitations in 
a normal metal have Floquet operator F* . Andreev refiections occurs at the opening 
to the superconducting reservoir, which is again represented by the matrix P. The 
symmetrized quantum Andreev map is finally constructed from the matrix product 

^ ^1/2/ F \_i/2 ^ / l-P^P -iP^P \ 



The excitation spectrum is obtained by diagonalization of JF, whose quasienergy spec- 
trum exhibits two gaps at E = and E = h/2TB. It can be shown that the excita- 
tion spectrum is identical to the solutions of the conventional determinantal equation 
det [1 -I- S{E)S*{E)] = 0, where the scattering matrix is given by Eq. (30). 
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